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Abstract. Confinement between two parallel surfaces is found, theoretically and 
experimentally, to drastically affect the hydrodynamic interaction between colloid 
particles, changing the sign of the coupling, its decay with distance and its 
concentration dependence. In particular, we show that three-body effects do not 
modify the coupling at large distances as would be expected from hydrodynamic 
screening. 



1. Introduction 

The dynamics of colloid suspensions and macromolecular solutions are governed by 
hydrodynamic interactions, i.e., correlations in the motions of particles mediated by 
flows in the host liquid 12] • I n various circumstances colloids are spatially confined 
by rigid boundaries as in, e.g., porous media, biological constrictions, nozzles, or 
microfluidic devices. Extensive studies have been devoted recently to colloid dynamics 
at the single-particle level, highlighting flow-mediated effects of the boundaries [H] 
Confined suspensions have been studied also by computer simulations [£] [Til]. 

Hydrodynamic interactions in an unconfined suspension jTTj decay with inter- 
particle distance r as 1/r. They are positive, i.e., particles drag one another in the 
same direction. The long range of the interaction leads to strong many-body effects, 
manifest, e.g., in an appreciable dependence of transport coefficients on particle volume 
fraction. A particularly important many-body effect is hydrodynamic screening [2] 
over length scales much larger than the typical inter-particle distance the suspension 
responds to a slow disturbance as if it were a homogeneous medium with merely an 
increased viscosity. This implies that at large distances many-body effects change the 
prefactor of the ~ 1/r pair interaction. 
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Our aim has been to investigate the effects of confinement on hydrodynamic 
interactions. Elaborating on our previous short publication |7], we consider here the 
case of confinement between two parallel plates, leading to a quasi-two-dimensional 
(Q2D) suspension (figure Q). In sectional we define the problem and the corresponding 
nomenclature. In section El we address the hydrodynamic interaction between two 
isolated particles, and in section 0] — the three-body correction at finite concentration. 
The experimental measurements are presented in section and in section |B] we discuss 
the results. 




Figure 1. (a) Optical microscope image of the experimental Q2D suspension at area 
fraction (f> — 0.338. (b) Schematic view of the system and its parameters. 



2. Definition of the problem 

The geometry considered in this work is depicted in figure ^6). Identical, spherical 
particles of radius a are suspended in a liquid of viscosity rj and temperature T, confined 
in a slab of width w between two planar solid surfaces. The x and y axes are taken 
parallel, and the z axis perpendicular, to the surfaces. The particles are assumed to 
behave as hard spheres with no additional equilibrium interaction. For simplicity we 
consider cases where particle motion is restricted to two dimensions, i.e., to a monolayer 
at the mid-plane, z = 0. We use the notation r(p, z) for three-dimensional position 
vectors, where p(x, y) is a two-dimensional position vector in the monolayer. The area 
fraction occupied by particles is denoted by 0. The Reynolds number is assumed very 
low. (It is of order 10~ 6 in the experimental system.) The hydrodynamics, therefore, 
are well described by viscous Stokes flows [O]. The confining boundaries are assumed 
impermeable and rigid, imposing no-slip boundary conditions on the flow. 

We characterize the pair hydrodynamic interaction between two particles by the 
coupling mobilities A LT (p) as functions of the inter-particle distance p. These are the 
off-diagonal terms in the mobility tensor of a particle pair, i.e., the proportionality 
coefficients relating the force acting on one particle with the change in velocity of 
the other. The two independent coefficients, Al and At, correspond, respectively, 
to the coupling along and transverse to the line connecting the pair. [In reference |Zj 
four coefficients were considered, A^ T , whose relation with the ones considered here is 
A LiT = (A+ T - A L>T )/2.] 
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We use dimensionless quantities throughout this paper, scaling all distances by 
the confinement length w and all mobilities by B a/w, where B Q = (G^r/a) -1 is the 
Stokes mobility of a single, unconfined sphere. In places where we refer to dimensional 
quantities for clarity, a tilde symbol is used, e.g., f = wr. 

3. Pair interaction at infinite dilution 

We first examine the effect of confining boundaries on the flow due to a local disturbance 
such as the one created by single-particle motion. The flow field can be found by solving 
the Stokes and continuity equations subject to the appropriate boundary conditions 
[TT] . This was done for a point force (stokeslet) in a Q2D geometry in reference |12j . 
However, by returning to the physical origin of the equations, one can gain a useful 
insight applicable to various geometries and particle sizes. The motion of a particle 
perturbs the local liquid momentum and displaces liquid mass. In an unconfined liquid 
the diffusion of the momentum disturbance, whose leading moment is a momentum 
monopole, creates a far flow decaying as 1/r. The mass term, whose leading moment 
must be a mass dipole (since mass is neither created nor lost), adds to the far flow a 
much smaller contribution ~ 1/r 3 . 

In confinement, however, the situation is very different. It is the momentum 
perturbation which, while diffusing away, is absorbed by the boundaries and thus gives 
a contribution to the far field which is exponentially small in r = r/w in both Q2D 
and Q1D geometries. The contribution from mass displacement propagates laterally, 
creating a far flow parallel to the boundaries. We conclude that the far flow field 
induced by single-particle motion in Q2D is that of a 2D mass dipole (source doublet). 
If the particle motion (mass dipole) is taken along the x direction, we thus have a far 
velocity field v(r) of the form 

v x (r) ~ f(z)(x 2 - y 2 )/p\ v y (r) ~ f(z)(2xy)/p\ v z (r) ~ 0. (1) 

The 2D dipolar field must be modulated by a transverse profile f(z) to make it vanish 
on the two confining surfaces. Note that the xy dependence of the far flow in Q2D, 
equation (0), differs from the far flow due to a local disturbance in a purely 2D liquid. 
The latter is governed again by (2D) momentum diffusion, leading to a logarithmic 
distance dependence. In ID only a mass monopole creates flow. Hence, in a Q1D 
geometry (confinement in a linear channel) both the momentum and mass contributions 
are exponentially small in r/w, leading to a short-ranged hydrodynamic interaction i6|. 

Exact calculations of the flow due to a point force (stokeslet), i.e., the Oseen tensors 
in the confined geometries, confirm the above conclusions ^21 EH] • Those exact results 
are relevant to colloid motion in the limit of very small particles, a <C w,r. The strength 
of the heuristic arguments presented above is their general validity — the conclusions 
concerning the far flow should hold regardless of particle size or details of the confining 
boundaries. These arguments also imply that, for a spherical particle, particle rotation 
does not contribute to the far flow under confinement, since it perturbs only the liquid 
momentum. 
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When the suspension is very dilute we may consider the interaction between two 
isolated particles. Suppose that the motion of particle 1 exerts on the liquid a unit force 
in the x direction. The disturbance will create a flow field which, in turn, will entrain 
particle 2 located at r. Thus, the coupling mobilities at large distances are obtained as 
Al = v x (px) and At = v x (py), where v(r) is given by equation (JTJ) and x, y denote unit 
vectors in the x and y directions, respectively. This yields 

0^0: A L (p>l) = A/p 2 , A T (p > 1) = -A/p 2 , (2) 

where A = f(z = 0) is a coefficient which depends only on a/w. (In the limit a/w <C 1 
one finds A = 9/16 L 12j.) As found in equation (J2J), the pair hydrodynamic interaction 
in Q2D is very different from its unconfined counterpart. The decay with distance 
is faster, ~ 1/r 2 instead of ~ l/r, yet the interaction is still long-ranged [Hj. (Its 
decay is slower than near a single surface, where the interaction falls off as l/r 3 [3J.) 
The transverse coupling is negative, i.e., particles exert "anti-drag" on one another 
as they move perpendicular to their connecting line. (In the unconfined case one has 
Ax = Al/2.) These properties are direct consequences of the far flow field inferred 
above. The 1/p 2 decay is that of the flow due to a 2D mass dipole. The negative 
transverse coupling is a result of the circulation flows in the dipolar field [7j. 

4. Three-body correction 

As cf) is increased, the pair hydrodynamic interaction should become affected by the 
presence of other particles. If hydrodynamic screening were to set in, we would expect 
the interaction at distances much larger than the typical inter-particle distance to have 
the form of equation (j2j) yet with a modified, ^-dependent prefactor. However, this is 
not the case, as is demonstrated below both theoretically and experimentally. 

We begin again with particle 1, located at the origin and exerting a unit force in the 
x direction. This creates a flow of the form (|Tj) . which entrains particle 2, located at p, 
with velocity v(p). We now introduce particle 3 at p' . Particle 3 will obstruct the flow, 
thereby exerting an extra force f3 on the liquid. This force is related to variations of the 
flow v(p') over the volume now occupied by particle 3. For example, in an unconfined, 
isotropic system, according to Faxen's first law [1, f 3 ~ V 2 v(p'). In the anisotropic 
Q2D geometry this law is modified to 

f 3i = C' ij Viv J (p',0) + C! J d zz \ z=0 v J (p',z), i,j = x,y. (3) 

[This result is readily obtained by inspecting all the ways to produce a vector f 3 from 
derivatives of another vector v while using the various symmetries of the field (0).] In 
equation © C, C" are 2x2 coefficient tensors (which can be found by a more detailed 
calculation), and = d xx + d yy . Since the flow (JTJ) satisfies V^v = and <9 2 | 2=0 v ~ v, 
equation (jHJ) is simplified to 

hi = C i3 v 3 {p'). (4) 

This local relation between f 3 and v must be invariant to rotation in the (x, y) plane, 
which leads to C xx = C yy = Cl and C xy = —C yx = Cr- Since the problem is linear, the 
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effect of the extra force on the velocity at p is found by projecting the same flow field, 
equation (JTJ), from p' onto p, Sv x = Vi(p — p')f 3i (p'). We then get 

SV X (P, p') = C L [V x (p - p')V X (p') +Vy(p~ p')Vy(p')] + 



Equation © presents the correction to the x-velocity of the liquid at p given a fixed 
position p' of the third particle. We are interested, however, in the correction averaged 
over all possible positions p' . This requires the probability density p(p') of finding 
a particle a distance p' away from particle 1, which we assume here to be uniform, 
p = (f>/(7ia 2 ). We average according to (6v x ) = j d 2 p'p6v x (p, p') and find 



The terms proportional to Ct cancel, leaving us with only one coefficient, e = Cl/ (ira 2 ). 

For p' ^> p the integrand in equation © decreases as (p')~ 4 - Hence, the integrated 
contributions from far-away third particles should lead to a correction proportional to 
4>/p 2 which, as expected, would renormalise the prefactor of the leading 1/p 2 term. 
However, carefully carrying out the convolution in equation (jHJ), as shown in the 
Appendix, proves this conclusion wrong. The specific angular dependence of the flow 
(JTJ) leads to cancellation of terms, and the integral of equation (0) vanishes. Thus, the 
obstructions from numerous particles add up to nothing, and the long-range response 
of the suspension is identical to that of the particle-free liquid (to linear order in <p). 
Equation (J2J), therefore, remains unchanged for low but finite concentrations. Deviations 
of the flow from the far-field dipolar form and the introduction of a nonuniform particle 
distribution (i.e., static pair correlations) lead to non- vanishing three-body effects. 
These corrections, however, are short-ranged and thus do not change the above result. 
(The short-ranged three-body effects will be addressed in a separate publication.) 

5. Experimental results 

The experimental system consists of an aqueous suspension of monodisperse silica 
spheres (diameter 2a = 1.58 ± 0.04 pm, density 2.2 g/cm 3 , Duke Scientific), tightly 
confined between two parallel glass plates in a sealed thin cell (figured). The inter- 
plate separation is w = 1.76 ± 0.05 pm, i.e., slightly larger than the sphere diameter, 
2a/w ~ 0.90. Digital video microscopy and subsequent data analysis are used to locate 
the centres of the spheres in the field of view and then extract time-dependent two- 
dimensional trajectories. Details of the setup and measurement methods can be found 
elsewhere |15j . Measurements were made at four values of area fraction, = 0.254, 
0.338, 0.547, 0.619 ± 0.001. From equilibrium studies of this system [TH] we infer that, 
for the purpose of this study, the particles can be regarded as hard spheres. 

The correlated Brownian motion of a particle pair is characterized by a longitudinal 
and transverse coupling diffusion coefficients. Upon scaling by D a/w, where D = 
k^T j (6nr]a) is the Stokes-Einstein diffusion coefficient of a single unconfined sphere, the 



C T [v x (p - p')vy(p') -v y (p- p')v x {p')\. 



(5) 




(6) 



Hydrodynamic interaction in quasi-two-dimensional suspensions 



6 




P P 

Figure 2. Longitudinal (Al, circles) and transverse (At, squares) coupling diffusion 
coefficients as a function of inter-particle distance p. The coefficients are scaled by 
Doa/w and the distance by w. Area fractions are (f> — 0.254 (a), 0.338 (b), 0.547 (c), 
and 0.619 (d). Dashed lines are a fit to ±A/p 2 with the same value of A = 0.36 for all 
panels. 

two coefficients (thanks to the Einstein relation) become identical to the two mobilities 
Al and At- The two coupling mobilities are thus directly extracted from measured 
particle trajectories as Al(p) = (xi(t)xz(t))p/(2Dot) and A-r(p) = (yi(t)y2(t))p/(2D t), 
where Xi(t) and Ui{t) are the displacements of particle i of the pair during a time interval 
t along and transverse to their connecting line, respectively. The average () p is taken 
over all pairs whose mutual distance falls in a narrow range (±0.09 /im) around p = wp. 

The measured coupling diffusion coefficients are presented in figure El The leading 
behaviour for all (ft values fits well the predicted ±A/p 2 dependence of equation (J2J). 
The negative sign of At confirms the predicted "anti-drag" between particles located 
transverse to the direction of motion. The fact that the fitted value of A = 0.36 does 
not change with demonstrates the absence of hydrodynamic screening up to an area 
fraction of 0.619. (Larger area fractions could not be checked because the suspension 
began to crystallise [IE]-) 
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6. Discussion 

Flows due to local disturbances in a Q2D geometry are governed by mass propagation 
rather than momentum diffusion. We have demonstrated the strong effect that this 
has on the pair hydrodynamic interaction between colloid particles. The 1/r 2 decay of 
the coupling with distance is faster than in the unconfined case but slower than near 
a single surface. The transverse coupling becomes negative. Three-body effects do not 
renormalise the coupling coefficient, in contrast with the usual case of hydrodynamic 
screening. Since momentum diffusion does not contribute to the large-distance coupling, 
it is perhaps natural that the flow effects of numerous distant particles may not 
renormalise the viscosity, which is the transport coefficient associated with momentum 
diffusion. 

We have not treated the effect of particle motion perpendicular to the boundaries. 
Such fluctuations must exist in practice, yet our experimental results suggest that they 
have a minor effect on the long-distance behaviour. This is expected since the flows 
produced by transverse fluctuations decay exponentially with distance j!2j . 

We have restricted the discussion to far-field effects. At small inter-particle 
distances various features are observed in the hydrodynamic interactions (see Fig. EJ), 
which are related to the static particle pair correlation. We shall address these short- 
ranged effects in a forthcoming publication. 
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Appendix 

We wish to calculate the convolution integral appearing in equation ©: 



Calculating I(p) by direct integration is quite tricky for reasons similar to those 
encountered in electrostatics. Excluding two small areas ~ a 2 around the (integrable) 
singularities at p' — and p' = p, we divide the integration into three domains (<p and 
if' being the polar angles of p and p', respectively): (i) a < p' < p — a, < if' < 2n; (ii) 
p' > p + a, < if' < 2n; (iii) p — a < p' < p + a, if + a/p < ip' < if + 2tt — a/ p. The 
contribution from the inner domain (i) vanishes upon angular integration for any a. The 



J(p) = J d 2 p'[v x (p')v x (p - p) + Vy(p')v y {p - p% 

where (omitting the prefactor) 

v x (p) = (x 2 - y 2 )/p\ v y {p) = 2xy/p\ 



(7) 



(8) 
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outer domain (ii) contributes ir/ p 2 +0(a/ p 3 ). However, integration over the intermediate 
narrow domain (iii) just cancels the outer contribution, yielding —Tr/p 2 + 0(a/p 3 ). Thus, 
in the limit a — > the result is / = 0. A finite value of a leads to a correction of 0(a/p 3 ), 
which is negligible compared to the leading 1/p 2 coupling at large enough distances. 

The direct integration pitfalls can be circumvented altogether by switching to 
Fourier space, jF[/(p)] = /(q) = J d 2 pe~ tcl ' p f(p). The transformed velocity field is 



Hence, / = for any p > 0. 
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v x (q) = -u{q 2 x - q 2 y )/q 2 , v v (q) = - 
Using equation (jHJ) in the convolution readily gives 
I(p) = n 2 J^- 1 [l] = n 2 5(p). 



n{2q x qy)/q 2 - 



(10) 
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